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We revisit the electromagnetic heat transfer between a metallic nanoparticle and a metallic semi- 
infinite substrate, commonly studied using the electric dipole approximation. For infrared and 
microwave frequencies, we find that the magnetic polarizability of the particle is larger than the 
electric one. We also find that the local density of states in the near field is dominated by the 
magnetic contribution. As a consequence, the power absorbed by the particle in the near field is due 
to dissipation by fluctuating eddy currents. These results show that a number of near-field effects 
involving metallic particles should be affected by the fluctuating magnetic fields. 



I. INTRODUCTION 



tromagnetic waves. A possible application to design 
near-field energy converters has been studied [33] . 



A lot of attention has recently been devoted to 
the interaction between nano-objects like atoms, 
nanoparticles or AFM tips, and surfaces, which is 
mediated by fluctuating thermal fields. A great 
variety of_ phenomena like Casimir-Polder forces 



[3, H U El, friction Forces 
field heat transfer 
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or near- 
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by the associated stochastic thermal currents. A 
common assumption is that the electric dipole ap- 
proximation can be used to model the nano-object 
[H El M, HI, [H, H, H, Here, we revisit 

the heat transfer between a surface and a metallic 
nanoparticle. We find that the leading mechanism is 
near-field induction heating, due to Joule dissipation 
of eddy currents in the particle. The large currents 
are produced by time-dependent infrared magnetic 
fields that dominate the energy density near a metal- 
lic surface. We find a different distance dependence 
of the flux as compared with the case of polar materi- 
als. Our work may find applications on local heating 
for data storage |27[ , and lithography [2(| HI] . 



The electric dipole moment of a sphere with radius 
R and dielectric constant e r is generally assumed 
to give the leading contribution [II El [H HH, I22I 
HEm, HH because it varies like (R/X) 3 whereas the 
next term in the Mie expansion varies as (R/X) 5 (X 
is the wavelength in vacuum) [34j . In this work, 
we will show that the interaction between the mag- 
netic dipole and the large magnetic fields in the near 
field may give the dominant contribution to the heat 
transfer. 

In the next section, we compare the absorption cross 
section of the electric and magnetic dipole moment. 
The third section is devoted to the analysis of the 
electric and magnetic energy density in the near 
field of a metal-vacuum interface. The final section 
analyses the heat transfer and discusses the physical 
mechanism. 



All the phenomena previously cited have to be de- 
scribed in the framework of fluctuational electro- 
dynamics introduced by Rytov 29]. It is known 
that the radiative heat flux between two bodies 
[H GH, El El El EE EE El El El can be dramat- 
ically enhanced when their separation distance be- 
comes smaller than 10/zm. It was found that evanes- 
cent waves yield the leading contribution to the heat 
flux. Experiments have been reported demonstrat- 
ing these effects [H HH, HH . It has also been pre- 
dicted that this heat transfer could have a very nar- 
row energy spectrum [lfj, [T7|, [2l| due to surface elec- 



II. ABSORPTION BY A METALLIC 
NANOPARTICLE 



Let us first compute the power absorbed by a small 
metallic particle. In what follows, we will use an 
isotropic, homogeneous and local form of the com- 
plex dielectric constant. To lowest order in R/X [34]], 
the particle can be described by its electric dipolar 
moment p. We define a complex polarizability ole 



p = a E e E, 



(f) 



where eo is the dielectric permittivity in vacuum and 
E is the external electric field. Another contribution 
is given by the magnetic dipolar moment fh charac- 
terized by its magnetic polarizability ajj 

fh = ajj H, (2) 

where H is the external magnetic field. Higher mul- 
tipoles can be neglected if \e r \ ^> 1 and R/X <C 
1. The contributions of the electric and magnetic 
dipoles to the power dissipated in the particle at a 
positive frequency u> are given by [HI, [3g] 

dH=^ImN e <|J? 2 |2> , (3) 



magnetic one because R/X <C 1. Yet, for values df 
e r such that |e, | ^> 1, as it is the case for metals at 
low frequencies, the magnetic dipole may provide the 
leading contribution. The physical reason is that the 
magnetic fields are continuous at an interface so that 
they can penetrate in the material. By contrast, the 
electric field in a spherical particle Ei nt is related to 
the external electric field by E int = [3/ (e r + 2)]E ext . 
Surface charges induced at the interface prevent the 
electric field to penetrate efficiently in the metallic 
particle. This screening effect takes place on a length 
scale given by the Thomas-Fermi length. It does not 
depend on the skin depth. 



P a t,H=w2Zm(a H ) Mo 



< \H\ 2 > 



(4) 



where /j,q is the magnetic permittivity in vacuum. 
ctE and an can be found in ref. I37f : 



a E = 47T R 



aH = TE R [— 



3 6 r 1 



(er - 1), 



(5) 



(6) 



where e r is the relative dielectric permittivity. 

Here, we do not take into account the diamagnetism 
of the material. Instead, the magnetic dipole mo- 
ment is due to eddy currents in the particle. The po- 
larisabilities arc calculated assuming that R is much 
smaller than the skin depth 5. A different form [35| 
can be derived when dealing with particles such that 
i<R<A. In what follows, we should keep in mind 
that the dipole model is a fair approximation pro- 
vided that the distance d between the center of the 
particle and a surface is much larger than R. Note 
that we have used for simplicity the extinction cross 
section of the elecric dipole. The exact form of the 
absorption cross section is discussed in ref. [H] . The 
difference for a metallic nanoparticle is negligible. 

As seen from Eqs (3,4), the absorption is the product 
of two terms, the imaginary part of the polarizabil- 
ity and the local density of energy. We shall show 
that for metallic nanoparticles at low frequencies, 
both terms are larger for the magnetic contribution. 
Let us first analyse the role of the polarizability. It 
appears from Eqs (3-6) that for values of the di- 
electric constant on the order of unity, the electric 
dipole contribution to losses is much larger than the 
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FIG. 1: Imaginary parts of the electric and magnetic 
polarisabilities of a gold sphere (lo p = 1.71 10 16 s -1 , vq — 
4.05 10 13 s" 1 , dp = 1.2 lOVs" 1 and A = 1). For u < v, 

Im(a H ) ^and for u> p , Im(a E ) ~^2L^l u 



We consider a non-magnetic metallic particle char- 
acterized by a Drude model e r = 1 — uj 2 / (uj 2 + i ujv) 
where uj p is the plasma frequency and v is the 
damping coefficient. To account for the confine- 
ment effects, the bulk dielectric constant e r [3§| 
is corrected by modifying the damping constant 
v = vq + A vp/R, where vq is the bulk damping 
coefficient, vf the Fermi velocity and A a sample- 
dependent coefficient. Figure 1 shows Im(etE) and 
Im{ptH) as a function of circular frequency for two 
gold spheres with radii R = 5 nm and R = 10 nm. 
It is seen that the electric polarizability is larger 
than the magnetic polarizability at optical frequen- 
cies. As explained before, this is no longer the case 
at low frequencies (typically smaller than v), where 
Im(a>H) is larger than 7m(as). 



III. LOCAL DENSITY OF ENERGY NEAR 
A METALLIC SURFACE 



To derive the energy absorbed by a particle in the 
vicinity of an interface, we need to consider the local 

densities of energy £0—^ — and po — —5 — ■ In a 
vacuum, both contributions are equal. The energy 
per unit volume U(z,ui) at a distance z from the 
interface increases dramatically in the near field due 
to the presence of evanescent waves as discussed in 
ref. [H, U(z,lu) is the product of the local 

density of states (LDOS) p(z, oj) by the mean energy 
of a mode 0(lu,T) = fku/[exp(huj/kBT) — 1], where 
2-kH is Planck constant, ks is Boltzmann constant 
and T the temperature of the substrate. The final 
expression for the evanescent part of the LDOS [H, 
|40| is the sum of the four following contributions : 
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where the subscripts E and M denote the electric 
and magnetic evanescent components, c is the light 
velocity in vacuum, p v (uS) — lu 2 /it 2 c 3 is the vacuum 
density of states, f(K,w) = 2(^-) 2 -l,r a = 

and r v = £l7o ~ e27 ° the Fresnel TE and TM reflection 
factors, and the complex number 7$ = 7^ + i 7^' is 
defined as the perpendicular part of the wave vector 

at a frequency u> : K 2 + 7q = £j 75- where i = in 
vacuum (eo = 1) and i = 1 in the metal. We have 
neglected non-local effects as we consider distances 
larger than the Thomas-Fermi screening length. 



On Figure 2, we plot the LDOS versus the fre- 
quency lu for d = 30 nm. The first conclusion is 
that the contribution due to the evanescent waves 
dominates. The second conclusion is that the s- 
polarized magnetic contribution is dominant for fre- 
quencies below lum — 2.4 10 14 s _1 which are relevant 
for heat transfer at 300 K. Indeed, in the expres- 
sion of the energy density U(z,lo), Q(lu,T) acts as a 
temperature-dependent frequency filter. At a given 
temperature, we define a cut-off frequency lum by 
ff M 6(lu,T)cLu/ J °° &(LU,T)dLu = 99/100. Frequen- 
cies much higher than lum are not relevant for heat 
transfer. We note that the p-polarized contribution 
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FIG. 2: Contributions of the evanescent waves to the 
local density of states (LDOS) at d = 30 nm of a gold- 
vacuum plane interface when using a bulk Drude dielec- 
tric constant. The magnetic and electric propagating 
contributions are also plotted (they are equal). 



associated with the surface plasmon polariton dom- 
inates at optical frequencies but does not contribute 
significantly in the infrared. 

The dominant contribution of magnetic energy can 
be understood by considering the analytical expres- 
sions of pf (Eq. 8) and p% (Eq. 9). Both ex- 
pressions are exactly symmetric, involving the same 
factor f(K, lu), and the imaginary part of the reflec- 
tion factor, respectively Im(r s ) and Im(r p ). The 
physical origin of the factor f(K,u) lies in a funda- 
mental difference of structure between propagating 
and evanescent waves. It is interesting to see why 
the magnetic energy dominates the electric energy. 
Indeed, in a vacuum, the electric and the magnetic 
energy are equal. This is no longer true for an s- 
polarized evanescent wave close to an interface. 



Let us denote the wave vector as follows 



fc = K e x - 7 e z 



(11) 



where K and 7 are the interface parallel and per- 
pendicular wave vectors and e x and e z unit vec- 
tors. The Helmholtz equation in a vacuum yields 
K 2 + 7 2 = lu 2 /c 2 . For an s-polarized field, the elec- 
tric field is given by 



E = (0,E,0), 



(12) 



and the magnetic field follows from the Maxwell- 

d 



Faraday equation in vacuum (curl E — — j^B) 



B=-(-7,0,K). 

LU 



(13) 



It follows that \B\ 2 = BB* = ^{\l\ 2 + K 2 ). For a 
propagating wave, this yields the well-known result 



\B\=c\E\, 



(14) 



whereas for an evanescent wave, 7 is purely imagi- 



nary so that |7| 



-7 . We get: 



c \E\ 




(15) 



For evanescent waves, K ^> fco so we find that the 
magnetic energy stored in an s-polarized evanes- 
cent wave is much larger than the electric energy. 
We have plotted in Figure [3] the function ^ = 

f{K,u>) = J i^ K ~ for different frequencies. It is 
seen that the density of energy, which is proportional 
to the density of states, is driven by the magnetic 
contribution. 
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FIG. 3: Ratio -jJj for s-polarized evanescent waves. 



Now, a similar reasoning can be done for p-polarized 
waves. In this case, the inverse ratio |-E|/c|I?| is also 
equal to <J f(K, lu) ~ \/2 K/ko which shows that 
the electric field dominates in this case. Near field is 
thus always dominated by an s-polarized evanescent 
magnetic field and a p-polarized evanescent electric 
field. The relative weight of both contributions is 
then given by the values of the imaginary part of the 
reflection factors. Since Im(r s ) is larger than Im(r p ) 
for a metal at low frequencies, the s-polarized contri- 
bution to the LDOS dominates as seen in Figure 2. 
As a consequence, the energy density is dominated 
by its s-polarized magnetic contribution. 



IV. HEAT TRANSFER BETWEEN AN 
INTERFACE AND A NANOPARTICLE 



We now combine the results obtained for the depen- 
dence of the polarisabilities and for the energy den- 
sity to derive the power absorbed by a small par- 
ticle as given by Eqs (3,4). Although heat trans- 
fer between a small particle and a substrate in the 
near field has only been calculated u sing the electric 
dipolar contribution 0, [H M, HI, HJ H, IU, 
the results shown in Figure 1 (large magnetic dipole 
moment) and in Figure 2 (large magnetic density of 
states) clearly indicates that the magnetic contribu- 
tion must be taken into account as suggested in ref. 
[32l |. Figure 4 shows the radiative power 



Prad — 



(j&M + i&H)^ (i6) 



dissipated by the substrate in the small particle 
(R = 5 nm). The key result observed in Figure 
4 is that the heat transfer is dominated by the s- 
polarized magnetic contribution. The magnetic con- 
tribution can be larger than the electric dipolar con- 
tribution by 3 orders of magnitude. The reason 
is that heat is dissipated essentially by eddy cur- 
rents. An important result is the dependence of the 
heat flux with distance. The magnetic LDOS varies 
asymptotically as \ jz and the electric LDOS as 1/z 3 
|36l . 140 ] . For gold, these behaviors are valid at very 
small distances (below 20 nm). Hence, there is no 
simple distance dependence for the absorbed power 
as seen in Figure 4. 
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FIG. 4: Radiative power dissipated in the gold parti- 
cle (radius R — 5 nm) by the semi-infinite planar gold 
substrate at 300 K, and the asymptotic behaviours 



The above analysis can be summarized by the fol- 
lowing scenario. Random currents flowing parallel 
to the interface can excite the s-polarized evanescent 
electromagnetic fields at infrared frequencies. As ex- 
plained above, the associated magnetic fields take 
large values in the near-field. They are continuous 
across a vacuum-metal interface so that they pene- 
trate efficiently in the nanoparticlc and can generate 
large eddy currents. These currents are dissipated 
through Joule effect. Thus, thermal heat transfer 
appears to be due to near-field induction heating. 
Radiative heat transfer between two parallel metallic 
surfaces can also be revisited with a similar scenario 




Heat is mainly dissipated by fluctuating eddy cur^ 
rents. The widely used electric dipole approximation 
is valid for dielectrics but breaks down for metals. As 
a consequence, the 1/z 3 dependence of the flux be- 
tween dielectrics is not valid for metals. A number 
of other effects due to thermal radiation (e.g. forces, 
friction) between metallic bodies are expected to be 
driven by their magnetic contribution, even if the 
media are non-magnetic. We note that the heat ex- 
changed by two metallic nanoparticles separated by 
a submicronic distance T^j should be driven as well 
by the interaction between their magnetic dipoles. 
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